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ABSTRACT 



Context. Heavy element abundances in extremely metal-poor stars provide strong constraints on the processes of forming these elements in 
the first stars. 

Aims. We attempt to determine precise abundances of strontium in a homogeneous sample of extremely metal-poor stars. 
Methods. The abundances of strontium in 54 very or extremely metal-poor stars, was redetermined by abandoning the local thermodynamic 
equilibrium (LTE) hypothesis, and fitting non-LTE (NUTE) profiles to the observed spectral lines. The corrected Sr abundances and previously 
^ obtained NLTE Ba abundances are compared to the predictions of several hypothetical formation processes for the lighter neutron-capture 
$_( elements. 

Results. Our NLTE abundances confirm the previously determined huge scatter of the strontium abundance in low metallicity stars. This 
' scatter is also found (and is even larger) at very low metallicities (i. e. early in the chemical evolution). The Sr abundance in the extremely 
, metal-poor (EMP) stars is compatible with the main r-process involved in other processes (or by variations of the r-process), as we briefly 
I ■ discuss. 

> 
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■rj" : 1. Introduction massive stars being the primary early producers of light ele- 

' ments (called a elements), this production declining at later 

^ Until recendy, our knowledge of the properties of the first stel- pj^^^^^ ^j^^ chemical evolution of the Galaxy, as (roughly) 

^ , lar generations of die Milky Way was quite restricted. Large ^^j. ^■^^ r-process elements 
J> ' telescopes and efficient detectors have made it possible to de- 

■r^ ■ termine the abundances of many chemical elements of very ^owcwc, developments in nucleosyndiesis models (e.g. 

^ primitive (very metal-poor) stars, constraining the theories of ™^Pold et al. M Fisher et al. M and Roberts et al. 

; nucleosynthesis in the first stars. 1™ show that the classical neutrmo-driven wind of core- 



The analysis of extremely metal-poor stars (EMP) provide 



collapse SNe is unable to produce the r-elements : they could 

. , , , , r , , instead be produced independently in other ways as proposed 

important information about the abundances of the elements , „, . , ^ . 

, , ^ J , , J- r bv Wanaio et al. (120101) . 

produced by the first stars, and about the corresponding forma- 

tion processes. The abundances of the heavy elements (Z > 30) The neutron-capture elements may also be produced by a 

are paiticulariy interesting. Several processes have been pro- "slow" neutron flux, but in the "main" s-process the produc- 

posed to explain flieir production, especiaUy two important tion takes place in intermediate-mass AGB stars, which have a 

neutron-capture processes: the rapid process (r-process) and the long evolution time, that is inconsistent with an eariy produc- 

slow process (s-process). Each of these processes is divided tion. Some early production could take place by means of the 

into two sub-processes ("main" and "weak"). "weak" s-process, in massive primitive stars although the lack 

Massive stars with a very short lifetime end their lives as °^ seed nuclei places this idea in doubt, 
core-collapse SN II, forming neutron stars, and it had been as- In the Solar System, the global abundances of the heavy el- 

sumed that the primitive massive stars provide a promising site ements are the result of a late production by the main s-process 

for the early production of the neutron-capture elements by the in addition to an early production by the r-process, and it has 

r-process (see e.g. Cowan & Thielemann ,2004, , Farouqi et al. been possible to distinguish the fraction produced by each pro- 

I2009I I. The hypothesis was thought to be consistent with these cess. 
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Cowan & Sneden (I2006I I demonstrated that the abundances 
of heavy (56 < Z < 70) elements, display nearly identical 
abundance patterns in the r-fraction of the abundances of the 
Solar System, and in 6 r-rich EMP stars (some more r-rich 
EMP stars have been added recently to the list). As a result, 
the main r-process has been proposed to be universal. With 
larger samples of analyzed stars, the whole picture appeared 
more complex. It is therefore quite important to determine the 
precise abundances of heavy elements in EMP stars. Since the 
electron concentration in the atmospheric plasma of these stars 
is substantially lower, significant deviations from LTE in these 
atmospheres are expected. Thus, reliable abundances of these 
stars can be obtained only by means of a NLTE analysis. As 
noted by Mashonkina & Gehren ( I2001I I. the NLTE correction 
can be significant, and clearly differs from one star to another 
We therefore decided to reanalyze the sample of the EMP stars 
previously studied by Frangois et al. (12007) and Bonifacio et 
al. (12009) in the framework of the ESO large program "First 
stars, first nucleosynthesis". For this work, we determined for 
these stars the abundance of Sr based on NLTE computations 
of the line profiles. We combine these abundances with previ- 
ous NLTE determinations of Ba abundances (see Andrievsky et 
al.. |2009 ) to study the evolution of the [Sr/Ba] ratio in the early 
phases of the chemical evolution of the Galaxy. 

2. Sample of stars and their parameters 

The spectra of the stars in this sample were presented in detail 
in Cayrel et al. ([2004] l and Bonifacio et al. (120071 12009I I. About 
fifty stars (35 giants and 18 tumoff stars), most of them with 
[Fe/H] < -3, were observed with the high-resolution spectro- 
graph UVES (Dekker et al., BOOOl l fed by the ESO-VLT The 
resolving power in the blue is R ^ 45000. The spectra were re- 
duced using the UVES reduction tool within MIDAS (Ballester 
et al.. 120001) . The signal-to-noise ratio (S/N) in the region of the 
blue strontium lines is typically ~ 120 per pixel with there be- 
ing an average of 5 pixels per resolution element. The sample 
includes a few well-known typical stars for comparison. 

The fundamental parameters (Tefj, log g, [Fe/H]) of the 
stars were derived by Cayrel et al. (12004) for the giants and 
Bonifacio et al. i2007) for the turn-off stars. The temperature 
of the giants is deduced from the colour by adopting the cal- 
ibration of Alonso et al. d 19961 1. and the temperature of the 
turn-off stars is deduced from the wings of the Ha line. The 
gravities are derived from the ionisation equilibrium of iron (in 
LTE approximation), and could be affected by NLTE effects. 
These parameters are repeated in Table [1] 

Briefly, MARCS model atmospheres have been used, but, 
for convenience, some parts of the computations have used 
ATLAS9 (Kurucz fT993l |2005|) model atmospheres wifli the 
overshooting option. It has been shown (Castelli et al. I1997I I 
that for EMP stars, non-overshooting models are more appro- 
priate. Non-overshooting Kurucz models have been shown to 
provide abundances very similar (within 0.05 dex) to those of 
the MARCS models used by Cayrel et al. (12004) and Bonifacio 
et al. (120091) . 

Table [U provides the adopted parameters and the Sr abun- 
dances. In addition we present the Ba abundances determined 



in the same way in Andrievsky et al. (120091 ) but using Kurucz 
models with the overshooting option (Kurucz Fl 9931 120()5T ). As 
a consequence, these previous NLTE Ba abundances are shifted 
(Table [TJ by -0.03 dex for the turn-off stars, and -0.05 dex for 
the giants. 

Three new giant stars were added to the original sample, 
their parameters and NLTE abundances (including barium) be- 
ing determined in the same way; they are listed at the end of 
Table [1] In the last column, we indicate the number of lines of 
Sr 11 measured. 

Two carbon-rich stars (CEMPs) are included in the sam- 
ple (they are designated by an asterisk in the first column of 
Table[B: CS 22949-037 (Sr-rich) and CS 22892-052 (a useful 
typical example of an extremely r-rich star). Among the turn- 
off stars, CS 29527-015 is a spectroscopic binary. 

3. Atomic model and calculations 

3.1. The strontium atomic model 

A NLTE analysis of the Srn spectrum has been the subject 
of only a few papers. A first application of the NLTE cal- 
culations for Ap-Bp stars was made by Borsenberger et al. 
( 119811 ). who used a six-level (plus continuum) atomic model 
of Srn. Belyakova & Mashonkina (1997) considered the sta- 
tistical equilibrium of Sr ii in atmospheres of F-G stars of dif- 
ferent luminosity classes with metallicity ranging from zero to 
-3. They used a model atom consisting of 40 levels of Srn 
and the ground level of the Srin ion, and this model was ap- 
plied to derive NLTE corrections. After this, Mashonkina & 
Gehren (120011 ) determined NLTE abundances in 63 cool stars 
with metallicity ranging from +0.25 to -2.20. Mashonkina et 
al. (120081 ) published results for the NLTE abundances of some 
heavy elements (including strontium) in four very metal-poor 
(VMP) stars with metallicity in the range from -2.1 to -2.7. 

Our strontium atomic model is quite similar to that first pro- 
posed by Belyakova & Mashonkina (119971) . The Grotrian dia- 
gram is given in Fig.[T] The model consists of 44 levels of Srn 
with n < 13 and / < 6, as well as the ground level of Srni (no 
excited levels of this ion were included, since they have very 
high excitation potentials). For the two important terms 4d^D 
and 5p^P °, the fine structure was taken into account. For the 
other levels, it was ignored. 

Since the ionization potential of the neutral strontium is 
only 5.7 eV, even in the atmospheres of cool stars it exists al- 
most only in the form of Sr ii. Therefore, the 24 levels of Sr i 
were included in the model only for the particle number con- 
servation. 

3.2. The calculations 

For all considered levels, the excitation energies were taken 
from Lindgard & Nielsen (119771 ). Oscillators strengths for 
72 bound-bound transitions were selected from three papers: 
Wiese & Martin (1980), Warner (1968), and Goldsmith & 
Boxman {19SI). For the remaining 258 transitions, we used 
oscillator strengths from Lingard & Nielsen (11977) . 
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Table 1. Parameters of program stars and strontium abundance. The Ba abundance has been corrected for overshooting. 



star 


Teff, K 


log g 


V„kms-' 


[Fe/H] 


6(Sr) 


[Sr/H] 


[Sr/Fe] 


[Ba/H] 


Rem 


tumofF stars 




















BS 16023-046 


6360 


4.5 


1.4 


-2.97 


-0.08 


-3.00 


-0.03 


- 


2 


BS 16968-061 


6040 


3.8 


1.5 


-3.05 


-0.46 


-3.38 


-0.33 


- 


2 


BS 17570-063 


6240 


4.8 


0.5 


-2.92 


+0.10 


-2.82 


+0.10 


-2.97 


2 


CS 22177-009 


6260 


4.5 


1.2 


-3.10 


-0.09 


-3.01 


+0.09 


- 


2 


CS 22888-031 


6150 


5.0 


0.5 


-3.30 


-0.13 


-3.05 


+0.25 


- 


2 


CS 22948-093 


6360 


4.3 


1.2 


-3.30 


-0.20 


-3.12 


+0.18 


-2.97 


2 


CS 22965-054 


6090 


3.8 


1.4 


-3.04 


+0.20 


-2.72 


+0.32 




2 


CS 22966-011 


6200 


4.8 


1.1 


-3.07 


+0.07 


-2.85 


+0.22 


-2.87 


2 


CS 29499-060 


6320 


4.0 


1.5 


-2.70 


-0.15 


-3.07 


-0.37 




2 


CS 29506-007 


6270 


4.0 


1.7 


-2.91 


+0.31 


-2.61 


+0.30 


-2.37 


2 


CS 29506-090 


6300 


4.3 


1.4 


-2.83 


+0.54 


-2.38 


+0.45 


-2.77 


2 


CS 29518-043 


6430 


4.3 


1.3 


-3.20 


+0.00 


-2.92 


+0.28 




2 


CS 29527-015 


6240 


4.0 


1.6 


-3.55 


-0.06 


-2.98 


+0.57 




2 


CS 30301-024 


6330 


4.0 


1.6 


-2.75 


+0.27 


-2.65 


+0.10 


-2.67 


2 


CS 30339-069 


6240 


4.0 


1.3 


-3.08 


-0.11 


-3.03 


+0.05 




2 


CS 31061-032 


6410 


4.3 


1.4 


-2.58 


+0.55 


-2.37 


+0.21 




2 


giants 


HD 2796 


4950 


1.5 


2.1 


-2.47 


+0.57 


-2.35 


+0.12 


-2.68 


4 


HD 122563 


4600 


1.1 


2.0 


-2.82 


+0.10 


-2.82 


+0.00 


-3.67 


5 


HD 186478 


4700 


1.3 


2.0 


-2.59 


+0.57 


-2.35 


+0.24 


-2.67 


4 


BD +17°3248 


5250 


1.4 


1.5 


-2.07 


+ 1.17 


-1.75 


+0.32 


-1.74 


2 


BD -18°5550 


4750 


1.4 


1.8 


-3.06 


-0.93 


-3.85 


-0.79 


-3.67 


3 


CD -38°245 


4800 


1.5 


2.2 


^.19 


-1.75 


-A.61 


-0.48 


^.72 


2 


BS 16467-062 


5200 


2.5 


1.6 


-3.77 


-2.30 


-5.22 


-1.45 


- 


2 


BS 16477-003 


4900 


1.7 


1.8 


-3.36 


-0.40 


-3.32 


+0.04 


-3.62 


2 


BS 17569-049 


4700 


1.2 


1.9 


-2.88 


+0.45 


-2.47 


+0.41 


-2.72 


4 


CS 22169-035 


4700 


1.2 


2.2 


-3.04 


-0.33 


-3.25 


-0.21 


^.17 


2 


CS 22172-002 


4800 


1.3 


2.2 


-3.86 


-2.03 


^.95 


-1.09 


^.82 


2 


CS 22186-025 


4900 


1.5 


2.0 


-3.00 


-0.18 


-3.10 


-0.10 


-2.94 


2 


CS 22189-009 


4900 


1.7 


1.9 


-3.49 


-1.28 


^.20 


-0.71 


^.62 


2 


CS 22873-055 


4550 


0.7 


2.2 


-2.99 


-0.40 


-3.32 


-0.33 


-3.37 


2 


CS 22873-166 


4550 


0.9 


2.1 


-2.97 


+0.27 


-2.65 


+0.32 


-3.54 


3 


CS 22878-101 


4800 


1.3 


2.0 


-3.25 


-0.58 


-3.50 


-0.25 


-3.57 


2 


CS 22885-096 


5050 


2.6 


1.8 


-3.78 


-2.00 


-4.92 


-1.14 


^.72 


2 


CS 22891-209 


4700 


1.0 


2.1 


-3.29 


-0.23 


-3.15 


+0.14 


-3.65 


2 


CS 22892-052* 


4850 


1.6 


1.9 


-3.03 


+0.47 


-2.45 


+0.58 


-2.25 


4 


CS 22896-154 


5250 


2.7 


1.2 


-2.69 


+0.67 


-2.25 


+0.44 


-2.35 


3 


CS 22897-008 


4900 


1.7 


2.0 


-3.41 


-0.08 


-3.00 


+0.41 


^.37 


2 


CS 22948-066 


5100 


1.8 


2.0 


-3.14 


-0.68 


-3.60 


-0.46 


-3.97 


2 


CS 22949-037* 


4900 


1.5 


1.8 


-3.97 


-0.88 


-3.80 


+0.17 


^.47 


2 


CS 22952-015 


4800 


1.3 


2.1 


-3.43 


-1.33 


-4.25 


-0.82 


^.57 


2 


CS 22953-003 


5100 


2.3 


1.7 


-2.84 


+0.30 


-2.62 


+0.22 


-2.45 


2 


CS 22956-050 


4900 


1.7 


1.8 


-3.33 


-0.80 


-3.72 


-0.39 


-3.95 


2 


CS 22966-057 


5300 


2.2 


1.4 


-2.62 


+0.12 


-2.80 


-0.18 


-2.95 


2 


CS 22968-014 


4850 


1.7 


1.9 


-3.56 


-2.17 


-5.09 


-1.53 




2 


CS 29491-053 


4700 


1.3 


2.0 


-3.04 


-0.33 


-3.25 


-0.21 


-3.74 


2 


CS 29495-041 


4800 


1.5 


1.8 


-2.82 


-0.13 


-3.05 


-0.23 


-3.28 


2 


CS 29502-042 


5100 


2.5 


1.5 


-3.19 


-1.98 


^.90 


-1.71 


^.62 


2 


CS 29516-024 


4650 


1.2 


1.7 


-3.06 


-0.63 


-3.55 


-0.49 


-3.67 


2 


CS 29518-051 


5200 


2.6 


1.4 


-2.69 


+0.45 


-2.47 


+0.22 


-2.82 


2 


CS 30325-094 


4950 


2.0 


1.5 


-3.30 


-2.37 


-5.29 


-1.99 




2 


CS 31082-001 


4825 


1.5 


1.8 


-2.91 


+0.52 


-2.40 


+0.51 


-2.17 


4 


additional giants 




















CS 22891-200 


4700 


1.3 


2.0 


-3.65 


-1.75 


-4.67 


-1.02 


^.28 


2 


CS 22949-048 


4800 


1.5 


2.0 


-3.25 


-1.39 


^.31 


-1.06 


^.33 


2 


CS 22950-046 


4650 


1.0 


2.2 


-3.56 


-1.12 


^.04 


-0.48 


^.49 


2 
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Fig. 1. The Grotrian diagram for Sr ii atom. 

Radiative photoionization rates for the s, p, and d levels are 
based on the photoionization cross-sections calculated with the 
quantum defect method and the corresponding tables of Peach 
( 1196 7 ). For / and g levels, we used the hydrogen-like approxi- 
mation (Lang [T998] l. 

Collisional rates for the transitions between the ground 
level 5s^S and 5p^P 5d^D, and 6s^S levels were estimated 
with the help of the corresponding formula from Sobelman et 
al. (I1981I I. For the remaining levels, the excitation rates by elec- 
tron collisions were approximated using the formula of van 
Regemorter (.1962) . Collisional excitations for the forbidden 
transitions were taken into account using the semi-empirical 
formula of Allen (1973 ) with a factor of 1. Collisional ion- 
ization caused by electrons was described using the Drawin 
( |1961| l formula. Inelastic collisions with hydrogen atoms were 
taken into account with the formula of Steenbock & Holweger 
dl 9921 1 with a scaling coefficient 0.01. This coefficient was de- 
rived from the fitting of the Sr ii line profiles in the solar spec- 
trum (Kurucz et al. |1984| ). 

The NLTE strontium abundance was determined with a 
modified MULTI code (Carlsson |1986l l. Modifications are de- 
scribed in Korotin et al. ( |19991 l. Since we use Kurucz's atmo- 
sphere models calculated with ATLAS 9 (Kurucz ,1993., 2005a . 
the necessary background opacities for MULTI are also taken 
from ATLAS 9. In the modified version, the mean intensities 
used to obtain the radiative photoionization rates, are calculated 
for the set of frequencies at each atmospheric layer, and then 
stored in a separate block, where they can be interpolated. This 
enables us to take account of the absorption in a great number 
of lines, especially in the UV region, which is important for the 
accuracy of the photoionization rate calculations. 

After the combined solution of the statistical equilibrium 
and radiative transfer equations, we derived the populations of 



Table 2. Parameters of the Srii lines. 





f 


log Trad 


log 7vw 


4077.7090 


7.063e-l 


8.130 


-7.70 


4161.7920 


1.945e-l 


8.110 


-7.60 


4215.5190 


3.303e-l 


8.100 


-7.70 


4305.4430 


1.811e-l 


8.080 


-7.60 


10036.6530 


1.219e-2 


7.340 


-7.63 


10327.3110 


7.393e-2 


7.320 


-7.63 


10914.8870 


5.754e-2 


7.270 


-7.63 



the 5s^S, Ad^D, 5p^P^, and 6s^S levels allowing us to re- 
produce the profiles of the lines of interest, provided that the 
broadening parameters were known. Those parameters were 
taken from the VALData-baseQ. 

In the visual part of the stellar spectra and near-IR range, 
there are only a few detectable lines of Sr ii, including the reso- 
nance lines 4077 A, 4215 A, and the subordinate lines 10036 A, 
10327 A, and 10914 A. We note that resonance lines are signif- 
icantly blended with iron, chromium, as well as some molec- 
ular lines. The IR lines are not blended. The 10914 A line is 
unavailable in the UVES spectra, but the two other IR lines 
10036 A and 10327 A are, even though in some BMP stars 
they are both detected at the noise level. 

To compare the NLTE strontium line profiles to the ob- 
served spectrum in the visual part of the spectrum one needs 
to use a combination of the NLTE and LTE synthetic spec- 
trum since the resonance Sr ii lines are to some extent blended 
with other metallic lines. This was achieved using the updated 
code SYNTHV (Tsymbal, 119961 ) which is designed for syn- 
thetic spectrum calculations in LTE. With this program, we cal- 
culated a synthetic spectrum for particular wavelength ranges 
including the Sr ii lines of interest taking into account all the 
lines from each region listed in the VALData-base. For the 
strontium lines, the corresponding /j-factors (factors of devi- 
ation from LTE level populations) that had been calculated in 
MULTI were included in SYNTHV, where they were used in 
calculating the strontium line source function (in particular, the 
first plot in Fig. |2]gives an example of this procedure applica- 
tion). 

3.3. Test calculations 

All the above listed input data and our strontium atomic model 
were checked using several test calculations. For this, we used 
high-resolution spectra of some stars with well known param- 
eters. 

First of all, we calculated profiles for the six Sr ii lines in the 
solar spectrum (visual and infrared range, Kurucz's et al. 1984 
solar flux spectrum) using the Kurucz (1996) solar atmosphere 
model with microturbulence of 1 km s The strontium lines 
parameters are given in Table |2] 

The Srii lines in the solar spectrum are closely repro- 
duced by our calculations (see Fig. |2]i with strontium abun- 
dance (Sr/H) = 2.92 (as recommended by Asplund et al. 12005 1 ). 

' http ://ams . astro. uni vie . ac . at/vald/ 
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10915.0 UA) 



Fig. 2. Profile fitting of the Srn lines in the solar spectrum. 
Observed spectrum - dots, NLTE profiles - smooth line, LTE 
profiles - dotted line. 

Table 3. Fundamental parameters of the reference stars and 
their NLTE Sr abundance. 



Star 


Teff, K 


log 8 


Vt, km S-' 


[Fe/H] 


(Sr/H) 


Arc turns 


4300 


1.5 


1.7 


-0.33 


2.26 


Procyon 


6500 


4.0 


1.8 


+0.00 


3.12 


HD 1581 


6000 


3.9 


1.0 


-0.22 


2.69 


HD122563 


4600 


1.1 


2.0 


-2.82 


0.10 



10 : 



1 - 



: -..Sr III 


■ 




SUN 










5s ^ 




:4d 









-4 -3 

log! 








4215.2 4215.5 4215.8 10036.1 10036.7 10037.2 10326.5 10327.0 10327.5)i(A) 



Fig. 3. The same as Fig. |2]but for HD 122563 spectrum. No 
LTE profiles are shown. 



For the sake of completeness, we also provide pure LTE pro- 
files for all investigated lines. 

We also compared the calculated Srii line profiles with 
those observed in the following well studied stars: Arcturus, 
Procyon, HD 1581, and the metal-deficient star HD 122563 
(which is also in our list of program stars). For these stars, 
spectra are available from the UVES archive of the Paranal 
Observatory Project (Bagnulo et al. l2003l l. These spectra were 
obtained with R = 80000 and have S/N of more than 300. Since 
the third Sr II IR line 10914A is unavailable in the UVES spec- 
tra (but is important for the elaborated strontium atomic model 
check), we used the Keck IR spectrum of HD 122563 previ- 
ously analyzed in Smith et al. (.2004) . Adopted parameters of 
these stars are given in Table |3] The estimated NLTE strontium 
abundance of these stars is given in the last column of this table. 

Figure [3] shows that the Srn lines in HD 122563 that be- 
long to the different multiplets are well reproduced with a sin- 
gle strontium abundance. 

Very good agreement is obtained between the abundances 
derived from the lines in the infrared and visible domains for 
the stars HR 1581 and HD 122563. For HD 122563, the abun- 
dance of Sr (and Ba) is also in very good agreement with the 
NLTE results of Mashonkina et al. (l2008l l. 



Fig. 4. The distribution of the ^-factors in the solar atmosphere. 



4. The behaviour of the Mactors 

Figure |4] shows the behaviour of the /j-factors for the solar at- 
mosphere. As one can see, ^-factors of the 5p level are smaller 
than those for 5 s level over the whole atmosphere. At the same 
time, one can note that since Srii is a dominating ionization 
stage of strontium at these temperatures, the population of the 
ground 5 s level does not deviate significantly from an equi- 
librium population. Thus, the resonance lines are stronger, but 
only slightly because in the range of their formation the b- 
factors of the 5i level are close to unity. 

For the solar metallicity, the lines 4077 A and 4215 A are 
only slightly affected by the NLTE effects (abundance correc- 
tion is about 0.01), while for all three IR lines, corrections are 
rather significant (about 0.3). 

For the giant stars with metal-deficient atmospheres 
([Fe/H= -3 and T^ff = 5000 K), the NLTE correction is 
somewhat different (Fig.|5]l. At low metallicity, all the lines are 
formed, on average, in the deeper atmosphere layers. The lower 
electron concentration leads to significant deviations from LTE 
in the level populations. The ground level of Srii is underpop- 
ulated at log T - -1, while excited levels are overpopulated. 
The pumping mechanism increases the population of 5p level, 
and consequently the populations of the higher levels become 
overpopulated too (since 4d level is more tightly connected to 
5p than to 5s, it is also overpopulated). 

For the infrared lines, the /j-factor of the higher level is 
smaller than that of the lower level and the Zj-factor of the lower 
level is larger than unity. Altogether, this means that IR lines 
tend to be stronger The corresponding NLTE corrections are 
about -0.5 dex and even larger depending on the metallicity of 
the atmosphere and other parameters. 

The behaviour of the resonance lines were first described in 
detail by Belyakova & Mashonkina d 19971 1, and later also con- 
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Fig. 5. The same as Fig.|4]but for the metal-poor atmosphere. 



sidered by Short & Hauschildt (I2006I I. Here we mention only 
two cases, which are graphically presented in Fig. |6] At the 
temperature of 6500 K and [Fe/H] lower than -2, these lines 
becomes weaker than those for LTE, and this effect progres- 
sively increases with decreasing model metallicity. This behav- 
ior reflects that at this temperature and metallicity Sr ii is no 
longer a completely dominant ionization stage (only up to 70% 
atoms persist in the form of a singly ionized strontium), and 
the ground level of Srii appears to be significantly underpop- 
ulated. For the lower temperatures, the ratio of Srm to Srii is 
about 0.01, and the ground level tends to have an equilibrium 
population. The change in the profile between NLTE and LTE 
is small, while the equivalent width is practically the same for 
both LTE and NLTE. This is because the line core is formed in 
the upper atmosphere layers where the /j-factor for 5 p level is 
much smaller than for 5i (NLTE profile is deeper than LTE), 
while the far wings are formed at the depth where the situation 
is opposite. This is a result of the difference in the flux redistri- 
bution between the line core and wings. Since in reality the far 
wings of the line are ignored, the direct use of the equivalent 
width value in NLTE analysis often produces not quite correct 
results. In the plot at the right in Fig. |6] we show the LTE and 
NLTE profiles of the line that comprise almost equal equiva- 
lent widths (as follows from the calculations), while the line 
core regions look different. Therefore, the correct means of de- 
riving NLTE strontium abundance should be individual profile 
fitting for each star using a complete NLTE computation. 

In Fig. |7] we show the dependence of NLTE-LTE abun- 
dance correction upon effective temperature and other parame- 
ters. We note that NLTE corrections of the strontium abundance 
depend on the effective temperature and gravity of the model, 
as well as the strontium abundance itself. 




Fig. 6. LTE (dotted line) versus NLTE (continuous line) profiles 
for the two temperature values. 
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parameters. 



5. NLTE strontium abundance determination in 
programme stars 

The NLTE abundances of strontium derived in our EMP stars 
are given in Table [1] and graphically displayed in Fig. 8. 
The NLTE barium abundances are also provided in the ta- 
ble (when available); they were computed with overshooting 
Kurucz models (Andrievsky et al. 2009) and are, to ensure ho- 
mogeneity with Sr abundances, shifted by -0.03 dex for the 
turn-off stars and -0.05 dex for the giants. 
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The Sr and Ba abundance results shown in the following 
diagrams are all NLTE abundances. 

This work confirms (with some improvements) the main 
trends previously found for the Sr abundances, determined as- 
suming LTE, especially those for the same sample of EMP 
stars: Fig. 8, is directly comparable to Fig. 12 (upper panel) 
of Bonifacio et al. ( I2009I I. A comparison of the figures (and 
of the tables) indicates that the NLTE corrections are moder- 
ate, but not negligible. The NLTE strontium abundances are, on 
average, about 0.2 dex higher than the LTE abundances, with 
variations from star to star 

It has long been known that the scatter in the Sr abundance 
is very large (e.g. Mc William. 1 1 998) . The scatter in the [Sr/Fe] 
ratios is slightly smaller when NLTE (rather than LTE) values 
are used, but even using the NLTE values, this scatter remains 
very large (Fig. 8). The scatter is also larger at low metallici- 
ties (e.g.: Mc William [19981 Frangois et al. 120071 Bonifacio et 
al. '2009' Lai et al. '2008): this trend seems to be confirmed in 
NLTE (Fig. 8), although we note that this conclusion depends 
on the position of the C-rich star CS 22949-037; since this is 
also an c-element-rich star and therefore does not have a very 
low metallicity (Aoki et al. ,2005) , its extremely low [Fe/H] 
value might not be representative. 

Despite the large scatter, the abundances of the EMP stars 
show that, as the iron abundance increases, the Sr abundance 
increases, in the mean, far more than Fe, the [Sr/Fe] ratio in- 
creasing rapidly at first (faster than a secondary process) and 
stabilizing at about the solar [Sr/Fe] ratio after some overshoot. 
The additional productions of Sr and Fe then remain at this 
constant ratio, up to metallicities of around about [Fe/H] - - 
1.5, and the scatter decreases (see Frangois et al., .2007j . 

The large enhancement of Sr, already detected for the LTE 
measurements (McWilliam 1998 Frangois et al.. 12007) in the 
(Ba-poor) stai- CS 22897-008 (diamond at left in Fig. 8) is 
confirmed by our NLTE determinations. This enhancement is 
larger than for any of those listed previously in "Sr-rich metal- 
poor stars" (Truran et al. 12002) . The star is intermediate be- 
tween the Sr-rich stars HD 122563 and the (even richer) star 
BS 16550-087 (Lai et al. l2008b . but it has a lower metallicity 
(and a lower Ba abundance) than the two Sr-rich stars. It could 
be, as well as HD 122563, considered as an example of a so- 
called "weak" r-process (Honda et al. 2007 1: its Eu abundance 
is too low to be measured but, by that very fact, could be com- 
patible with such a "weak" r-process. 

The low metallicity turn-off star CS 29527-015 (extreme 
left closed circle in Fig. 8) is Sr-rich, its Ba and Eu abundances 
being too low to be measured. 

5.1. Srvs. Ba 

Fig. 9 shows [Sr/Fe] versus [Ba/Fe]. Both Sr and Ba globally 
increase with Fe. Since the Sr/Fe ratio has a larger scatter at low 
metallicities, it also has more scatter at low values of [Ba/Fe]. 

In this figure, the thick line corresponds to the constant ra- 
tio [Sr/Ba] - - 0.31, which is the value chosen by Qian & 
Wasserburg (12008) to represent their component for the Sr- 
poor stars. This ratio is compatible with the production ratio 




-4 -3.5 -3 -2.5 -2 



[Fe/H] 

Fig. 8. [Sr/Fe] vs. [Fe/H] in our program stars: dwarfs: 
closed circles, giants: open circles, the C-rich stars CS 22949- 
037 and CS 22892-052: asterisks, the Ba-poor star CS 22897- 
008: diamond. Three newly analyzed giants are added. Typical 
error bars have been indicated. 



of Sr and Ba for the main r-process: it is near (but very slightly 
smaller than) the ratio given by Burris et al. (2000) for the r- 
process-only component of these elements in the Solar system, 
it is also slightly below the similar value provided for this com- 
ponent by Arlandini et al. (I1999I I. 

In Fig. 9, all the stars (apart from one, marginally) are sit- 
uated higher than the line suggesting that, for most stars, one 
(or more) additional production(s) of Sr occurred (or a varia- 
tion in the main r-process, or both). The departures from the 
main r-process are higher at low metallicities (and at low Ba 
abundances). 

In our sample, two r-rich stars (CS 31082-001, CS 22892- 
052) defining the main r-process (Cowan & Sneden. 12006) . are 
(as expected) near the solid line (and also BD -i-17°3248). In 
their survey of r-rich stars (r-II stars), Barklem et al. (12005) note 
that the r-II stars are centred on a metallicity of [Fe/H] = -2.81, 
and the two r-rich giants in our sample, have indeed such rela- 
tively high metallicities (and are consequently relatively Sr and 
Ba-rich, in the figure being near the upper corner). The Sr-rich 
turn-off star (closed circle at left in Fig. |8) has a significanly 
lower metallicity. 

In Fig. 9, it is seen that the upper envelope has a [Sr/Fe] ra- 
tio that is approximately constant, while [Ba/Fe] increases by 
an order of magnitude, suggesting that the production process 
of Sr is different (somehow independent?) of the process pro- 
ducing Ba. With increasing Fe, [Ba/Fe] increases thanks to the 
main r-process, while along the upper envelope [Sr/Fe] does 
not. This at least indicates that an early production of Sr oc- 
curred in the chemical evolution of the Galaxy. 

Fig. 10 shows the positions of the stars in the [Sr/Ba] vs. 
[Fe/H] plot, together with a plot of the theoretical LEPP model 
predictions (Travaglio et al. 12004) . The corresponding theoret- 
ical curve is not far from the upper envelope of the data points, 
and an inhomogeneous mixing of the products of such a LEPP 
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Fig. 9. [Sr/Fe] vs. [Ba/Fe]. The solar ratio [Sr/Ba] = is rep- 
resented by the dashed line. The value [Sr/Ba] - -0.31, cor- 
responding to the main r-process production, is represented 
by the full line. Closed circles: dwarfs, open circles: giants, 
crosses : three added EMP giants. The stars without Ba mea- 
surements are not represented. Symbols as in Fig. 8. In Fig. 9 
and in following figures, the error bars are no longer indicated 
for an easier visual perception of the abundance trends. 
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Fig. 10. [Sr/Ba] vs. [Fe/H].LEPP model prediction represented 
by the thick line, [Sr/Ba] - -0.31 : dashed line, linear regres- 
sion : dotted line. 

process with the products of the main r-process could explain 
the distribution of the stars in this figure. 

In our sample, the carbon-rich EMP star CS 22949-037 is 
classified as a "CEMP-no" star: this classification is explained 
e.g. by Sivarani et al. (2006). This carbon-rich star with the 
rather high ratio [Sr/Ba] ^ 0.6 could also be explained by the 
LEPP process (or by massive rotating stars e.g., Pignatari et 
al. 120081 Meynet et al. l2010l or other processes; see discussion 
below). 

5.2. Production of Fe, Sr, and Ba 

Fran9ois et al. (I2007I I demonstrated (their Fig. 15) that in ex- 
tremely metal-poor stars, the ratio [Sr/Ba] (as well as Y/Ba and 
Zr/Ba) is often higher than solar in the range -3.7 > [Ba/H] > 
-A.l. A correlation is found between [Sr/Ba] and [Ba/H] in the 
interval -5.5 <[Ba/H]< -4.5 (LTE values). At very low metal- 
licity, Frangois et al. (2007) find that, for the lower observed 
values of [Ba/H] (about -5 dex), the ratio [Sr/Ba] returns to the 



Fig. 11. [Sr/Ba] vs. [Ba/H]. [Sr/Ba]H = -0.31 - dashed line, 
linear regression - dotted line. 

solar values. Our NLTE homogeneous determinations qualita- 
tively confirm this behaviour, and enable one to robustly claim 
that the Sr abundances are generally higher than those predicted 
by the main r-process pattern. 

Fig. 11 shows that, with [Ba/H] increasing from -5.0 to - 
4.0, the upper envelope of the scatter diagram reaches the one 
to one slope when excluding the star CS 22897-008 (and even 
more when including this star). Along this envelope, the Sr 
abundance increases much more than the Ba abundance. In the 
other part of the figure, where [Ba/H] increases from -3.5 to - 
2.5, the upper envelope decreases with a slope around -0.7; the 
Sr abundance increases more gradually than Ba. This behaviour 
illustrates the well-known lag between the onsets of Sr and Ba 
production, again suggesting that the production processes of 
Sr and Ba are different. 

Many stars are relatively near the dashed line (correspond- 
ing to the main r-process), about half of the stars having Sr/Ba 
ratios smaller than three times the main r-process ratio. 

The three stars with the lowest Ba abundances (including 
CD -38°245) have Sr/Ba ratios near the main r-process ratio (in 
agreement with the behaviour of the most Ba deficient star of 
Lai et al. (2008 ). 

Fig. 12 shows the simultaneous enrichment of Fe, Sr, and 
Ba, the low metallicity stars being both Ba-poor and Sr-poor. 
Fig.[T2]is somewhat similar to Fig. 9 but provides an indication 
of the iron production. 

Summarizing, we can state that our NLTE Sr and Ba abun- 
dances discussed above appear to be consistent with the Qian 
& Wasserburg (2008 1 simple three-component model of chem- 
ical evolution. According to that model, the first phase of the 
chemical evolution of the Galaxy would be the production of 
elements by high-energy hypernovae: their yields would in- 
clude significant amount of low-A elements (from Na to Zn, 
including iron), but no heavier species. In reality. Fig. 12 shows 
that iron-poor stars have a very low content of Sr and Ba. The 
shortfall in Sr later decreases to zero when the low-mass and 
normal-mass SNe {H and L sources respectively) begin pro- 
ducing strontium (Fig. 8). Moreover, since Sr is produced by 
both the H and L sources, and Ba is produced predominantly 
by the lower-mass H source, it can explain the general decrease 
in the Sr/Ba scatter, which is clearly seen seen in Fig. 11. 
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-fi n l- H 

-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 [ga/H] 

Fig. 12. [Sr/H] vs. [Ba/H]. Stars with [Fe/H] < -3.5 - 
small circles, —3.5 < [Fe/H] < -3.0 -intermediate circles, 
[Fe/H]> -3.0 - large circles. 

5.3. Comparison to faint dwarf spheroidal galaxies 

We have much to learn about dSphs, but a full discussion is be- 
yond the scope of this paper. We present only two short com- 
parisons. 

Cohen & Huang (1201 Oi l find, in the (metal-poor) stars of 
the Ursa Minor dwarf spheroidal galaxy, the heavy elements 
produced follow the main r-process pattern: this is consistent 
with the short duration of the star-forming phase of this small 
galaxy, inferred from the CMDs (Orban et al. 12008) . 

Norris et al. (12010) identify CEMP-no stars (or similar) in 
dwarf galaxies, and note a similarity between their abundance 
patterns and those of the same stars in our Galaxy: the similar- 
ity is indeed striking, but strontium is systematically stronger 
in the stars of our Galaxy. The comparison is made for LTE 
abundances, but for one of them (CS 22949-037), the high Sr 
abundance is confirmed by our NLTE determinations. 

5.4. Sr formed by other processes? 

Many independent explanations of excess of Sr production rel- 
ative to the main r-process have been proposed, such as 

1) a quenched flux model (Tmran et al. l2002l l: 

2) particular models of hypernovae (e.g. special computa- 
tion of four models of hypernovae reproducing the high Sr/Ba 
ratios (and other abundances) in four stars of our sample, by 
Izutani et al. ( l2009l : 

3) the inhomogenous mixing of three types of nucle- 
osynthesis productions: a phenomenological "three-component 
model" ( Qian & Wasserburg|2008ll; 

4) a weak r-process (Honda et al. I2007I I or a variation in the 
r-process (Roederer et al. 12010) : 

5) contamination by main s-process (by an AGB compan- 
ion): metal-poor AGB do not generally provide a strong Sr ex- 
cess, still, monitoring a subsample of Sr-rich stars would be 
useful; 

6) a LEPP process (Travaglio et al. l2004l l. 

Since nitrogen is observed as a primary product in some 
EMP stars (Spite et al. 2005), rotating massive stars may pro- 
duce heavy elements, and in some cases strong Sr/Ba ratios 
(Pignatari et al. I2008I I. However, the correlation of Sr/Ba with 



the abundance of N is weak, suggesting that, if rotating stars 
are a source of Sr, it is not the unique additional source. The 
Sr-rich star CS 22949-037, classified as CEMP-no, might be 
produced by a LEPP, or by a massive rotating star (producing 
also carbon). 

6. Conclusion 

The homogeneous NLTE Sr abundances of the stars of this 
sample, provide precise information about the Sr abundances 
in the early evolution of the Galaxy. In addition to previously 
determined NLTE Ba abundances, these data provide some in- 
sight into the process(es) of Sr production. 

The methods of the NLTE determination of Sr abundances 
have been verified (solar spectrum, extension in infrared), 
the results being in satisfactory agreement with computations 
found in the literature (e. g. Belyakova & Mashonkina 1 9971 
Mashonkina et al. 120081) and comparisons with a few stellar 
NLTE abundances available in the literature show good agree- 
ment. 

Our study here has confirmed the behaviour of the Sr LTE 
abundances found in the literature for both metal-poor stars 
(e.g. by McWilUam. [19981 and Lai et al. [2008] l and our sam- 
ple (Frangois et al. l2007l and Bonifacio et al. l2009l l. These be- 
haviours are now established on a firmer basis. 

During the chemical evolution of the Galaxy, when the 
iron abundance increases, the mean Sr abundance increases, the 
mean [Sr/Fe] ratio increases, reaching (and then exceeding) the 
solar ratio, before decreasing slowly towards the solar ratio. 

The Sr/Ba ratio is always higher than the main r-process ra- 
tio. In addition to the large scatter in the Sr/Ba ratios, this sug- 
gests that either some production of Sr (at least partly indepen- 
dent of the production of Ba) is added to the main r-process, or 
that the r-process that produces Sr and Ba varies significantly. 

A number of additional processes of Sr enrichment have 
been proposed in the literature, although none have been clearly 
verified, the data even being consistent with several processes 
existing simultaneously, enhancing the difliculty of identifying 
these processes. 

In this complex situation, it is of the upmost importance to 
collect new accurate abundances for a larger sample of stars, 
and a larger number of heavy elements (first and second peak 
elements) in the metallicity range -4.5 < [Fe/H] < -2.5, to try 
to identify the processes (and their sites) responsible for pro- 
ducing these elements. 
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